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(— i We describe some of the recent progress in our understanding of the dynamics of strongly cor- 

Qh related gluons at high parton densities. Computations in the Color Glass Condensate effective 

Q_i field theory provide a good description of inclusive and semi-inclusive final states in DIS, p+p 

and p+A collisions at small x. In nucleus-nucleus collisions, they provide an ab initio description 
of entropy generation, decoherence, isotropization and the onset of hydrodynamic flow. The suc- 
i—\ cessful description of a wide range of phenomena from RHIC to LHC is outlined and possibilities 

for more stringent tests noted. 
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Why is it of fundamental importance to understand the properties of strongly correlated gluons 
at very high energies ? Simply, it is because it offers a promising route to understand empirically the 
fundamental structure of matter. Despite QCD being established as the fundamental theory of the 
strong interactions, we are still far from having clear pictures of hadrons and nuclei that transcend 
even limited contexts. For instance, the constituent quark model (and variants) gives quite a good 
picture of the spectrum of hadron resonances, but doesn't do as well with the internal spin dynamics 
of the proton, and has little of value to say about high energy scattering. The sophistication attained 
by lattice QCD is impressive to behold but there is scant hope that it too will have much to say about 
multiparticle production in the deep Minkowski region of the theory. 

At the heart of the matter (no pun intended), are the outstanding problems of confinement 
and chiral symmetry breaking (and their interplay) that continue to thwart our developing a useful 
unified picture. In particular, there is an unresolved tension between the QCD parton model and the 
constituent quark picture. The former is not much use for spectroscopy but provides the framework 
for spectacularly successful computations of very high transverse momentum processes with all the 
"unpleasant" stuff factorized into parton distributions. These are universal and therefore provide 
predictive power. Perturbative computations however only reliably describe a small fraction of high 
energy cross-sections with the bulk parametrized by "stringy" event generators. The latter, albeit 
extremely useful (in particular for searches beyond the standard model), at best provide limited 
insight into hadron structure. 

Genuine progress in understanding hadron structure can be achieved in the Regge-Gribov 
asymptotics of very high energies at fixed but large momentum resolution scales. In these asymp- 
totics, a key feature of perturbative QCD is the rapid bremsstrahlung growth of gluon distributions. 
At fixed transverse resolution Q 2 of a probe, at a given parton momentum fraction x, this growth 
generates a phase space occupancy in the hadron approaching the maximal parametric value of 
1 /cCs permitted by the stability of the theory. For each value of Q 2 , there is a corresponding x value 
where gluon distributions saturate for momenta below a saturation scale Qs{x) corresponding to 
maximal occupancies. This phenomenon called "gluon saturation" [1] appears to be a fundamental 
feature of the theory 1 . In addition, the dynamical generation of a large scale Q$ 3> Aqcd suggests 
that the very non-perturbative properties of the hadron and hadronic scattering can be understood 
in weak coupling because GtsiQs) <C 1 in Regge-Gribov asymptotics. Furthermore, since the phase 
space occupancies are large (l/(Xs(Qs) ^ 1)> the saturation regime is a classical regime of the 
theory [2]. 

Hadron and nuclear wavefunctions are many body configurations of (primarily) gluons in this 
high energy limit. Since we have a weak coupling theory, one can properly ask what the effective 
degrees of freedom are, what is the nature of their correlations, how precisely does the coupling 
run, and whether there is a universal fixed point common to hadrons and nuclei. Ultimately one 
would like to understand how this weak coupling yet non-perturbative dynamics at distance scales 
l/Qs is altered at larger distance scales by chiral symmetry and confinement. For a number of 
bulk observables, if Q$ is large enough, can one anticipate that weak coupling computations can be 
performed with minimal contributions from these soft dynamics ? 

'The phase space occupancy argument is specific to the light front parton model formulation of the theory. The 
arguments can be formulated as well for gauge invariant quantities such as field strengths squared or "dipole" scattering 
amplitudes in the same framework. 
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The features mentioned as well as the kinematics of the high energy limit can be formulated 
in a weak coupling effective field theory (EFT) that describes the physics of saturated gluons as 
a Color Glass Condensate (CGC) [3]. The degrees of freedom 2 in this EFT are dynamical gluon 
fields at small x that couple coherently to static classical color charges at large x. Wilsonian renor- 
malization group (RG) equations, derived from requiring that observables be independent of the 
separation in x between sources and fields lead to an infinite hierarchy of evolution equations in x 
for n-point Wilson line correlators averaged over dense color fields in the nucleus. For a physical 
observable defined by an average over all static color source configurations, 



J[Dp}W Y \p}0\p], (D 



one can define a gauge invariant weight functional ("density matrix") of color sources W Y [p] at 
rapidity Y = ln(x/^o) = hi^o A~)> where xq is the initial scale for small x evolution, satisfies the 
JIMWLK equation [5] dW Y [p]/dY = Jt?W Y [p]. From eq. (1) then, the energy evolution of the 
observable can be expressed as 

d ML = {j^e) y . ( 2) 

The structure of the JIMWLK Hamiltonian 3%? is such that Y is an object distinct from (&} Y ; 

one obtains an infinite hierarchy of evolution equations for operator expectation values ( ff) Y . Given 
appropriate initial conditions at large x, solutions of this Balitsky- JIMWLK hierarchy [6, 5] allow 
one to compute a wide range of multi -particle final states in deeply inelastic scattering (DIS) and 
hadronic collisions. A remarkable result is that the full dynamical content of this infinite hierarchy 
of correlation functions can be obtained explicitly numerically by Langevin techniques [8] as has 
been demonstrated in computations we shall describe shortly. 

The simplest example of high energy scattering observables are the inclusive DIS structure 
functions F2 and Fi that are proportional to the forward scattering amplitude of a qq "dipole" on a 
nucleus. The forward dipole amplitude (dipole cross section) is given by 



a dip .(x,r T )=2 J d 2 b r x^l-— TrV(b r + yJV T (br-y M, (3) 

where rj = xj — yr is the transverse size of the dipole, bj = (xr + yr )/2 is the impact parameter 
relative to the hadron and V(xt) = Pexp(ig f dz~ P yf )■ In general, the Balitsky-JIMWLK 
equation for the expectation value 5 = (5), with 5(xr — yr) = 77: TrV(x7)V t (yr), also depends on 
the correlator (S(xj — zt)S(zt — yr))- In the large _/V c limit and for large nuclei, this factorizes to 
the form 55; numerical simulations of the B -JIMWLK hierarchy demonstrate this approximation to 
be much better than 1 /N c 2 in all but a very small kinematic window [9, 10]. With this factorization, 
one obtains the closed form Balitsky- Kovchegov (BK) equation [6, 7] 



^5(x r -yr) = yjr^ ^, yz [5(x r -z r )5(z r -yr) -5(x T - yr)] • (4) 



where J^ yz = (xr — yr) 2 /(xr — zr) 2 (zr — yr) 2 - A significant amount of work is underway in 
performing next-to-leading log x extensions of the BK equation [ll]-for the state of the art, we 



2 There is an increased understanding now of the dictionary between these d.o.f and those in the Reggeon Field 
Theory framework [4], 
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refer the reader to ref. [12] and references therein. At the phenomenological level, a running 
coupling BK equation has been developed (rcBK) [14], which gives good agreement with the latest 
combined HI and ZEUS data sets [15], and does better than NNLO DGLAP fits in the small x, 
small Q 2 region of the data [16]. 

For less inclusive observables, one encounters the expectation value Q = (Q) of the "quadrupole" 
operator 

<2(yr,xr,U7-,v T ) = ±-Ti (V + (yr)V(x r )V t (u 7 -)V(v 7 -)) . (5) 

Unlike (55), it is not reducible to the product of dipoles even in the large N c and large A approxima- 
tions and is a novel universal correlator in high energy QCD [17, 18], interesting both from theoret- 
ical and phenomenological perspectives. Examples of where this quantity appears include small-* 
di-jet production in e+A DIS [18], quark-antiquark pair production in hadronic collisions [19] and 
near-side long-range rapidity correlations [20]. Another interesting quantity [18] is a six-point cor- 
relator 56 (QS) that is probed in forward di-hadron production d + A — > h\ \i2 X in deuteron-gold 
collisions at RHIC. The dominant underlying QCD process is the scattering of a large x\ valence 
quark from the deuteron off small xi_ partons in the nuclear target, with the emission of a gluon 
from the valence quark either before or after the collision. The RHIC experiments show that the 
away-side peak in the di-hadron correlations is significantly broadened for central collisions at 
forward rapidities [21] as predicted in the CGC framework [22] and confirmed in more detailed 
analyses [23]. 

These analyses however relied on simplified factorization assumptions that are not justified. 
The JIMWLK RG equations for Q [24] and 56 [25] have been derived and computed numeri- 
cally for particular configurations of these operators [10]. A noteworthy result of these numerical 
simulations is that the initial condition and the RG evolution of higher point correlators are well 
reproduced by expressions that are functions only of the dipole expectation value 5. This result 
indicates strongly that the B-JIMWLK hierarchy is closely approximated by a Gaussian effective 
theory for W[p] whose variance is proportional to 5: the energy evolution of the theory is con- 
trolled entirely by BK-evolution [13]. These results are illustrated in fig. 1; the left figure shows 
the spatial (dipole) correlator of Wilson lines for a particular configuration of color charges from 
the numerical simulations of JIMWLK. The right part of the figure shows a comparison of the nu- 
merical results to the Gaussian approximation. Excellent agreement is observed for the expectation 
value 56 of the operator appearing in di-hadron correlations, while the factorized expression (5) 3 
disagrees strongly with the numerical results. 

The framework we have outlined here in terms of "dipole" and "quadrupole" operators has 
been used recently in phenomenological computations to describe the recent d+Au RHIC data [26]. 
More detailed computations and predictions for p+Pb results at the LHC are desirable [27]. 

As this brief discussion suggests, the CGC approach to gluon saturation has been very fruitful 
with rapid ongoing developments in our understanding of the structure of matter at high energies. 
How does one best relate these developments to the outstanding issues raised at the outset ? There 
are many theoretical approaches feasible, but to my mind high luminosity studies of exclusive 
processes at a future Electron Ion Collider [28] offer the greatest promise of unraveling the spatial 
dependence of high energy evolution. This relates the weak coupling non-perturbative dynamics 
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Figure 1: Left: Correlation l/A^ c (V^(0,0)y(x,y)) between the center position (0,0) and the point (x 7 y) 
for a given rapidity Y. Colors illustrate the degree of fluctuations. From ref. [10]. Right: JIMWLK nu- 
merical simulation for a particular "line" configuration of the six point Wilson line correlator (§(,) versus 
1 — (S), compared to an analytic Gaussian mean field approximation (MFA) for six different rapidities. 
From ref. [13]. 



of gluon saturation to the intrinsically non-perturbative dynamics of chiral symmetry breaking and 
confinement, as exemplified in the high energy limit by the Froissart bound [29]. 

We shall now discuss some of the implications of the the weak coupling ideas of the CGC EFT 
for high energy hadron scattering, in particular for nucleus-nucleus collisions. In this framework, 
the collision can, to leading order, be described as a collision of "shock waves" of classical gluon 
fields [30]. The preequilibrium matter produced is called the Glasma [31] and its properties can 
be determined ab initio in the CGC EFT framework. The classical field configurations are highly 
unstable; small quantum fluctuations around the classical fields are unstable and grow nearly ex- 
ponentially [32]. This rapid growth can cause the highly anisotropic initial Glasma color field 
configurations to isotropize on extremely rapid intervals not significantly greater than \/Q s [33]. 
This picture of underlying expanding classical fields and rapidly growing quantum fluctuations, 
whose interactions with the classical field isotropize the fields on very short time scales, bears 
strong analogy to "pre-heating" inflationary dynamics in the very early universe that generates the 
"hot" era of thermalized matter [34]. In particular, in direct analogy to the inflationary framework, 
preequilibrium dynamics is driven by a spectrum of amplified initial quantum fluctuations [35, 48] 
potentially generating a significant amount of flow before thermalization. 

A master formula which incorporates the leading divergent contributions from quantum fluctu- 
ations, both before and after the collision, has been derived to leading accuracy in a weak coupling 
expansion of inclusive quantitities; for the components of the stress-energy tensor, one obtains [36] 

(rHLLx+Linst. = / i D Pi D P2] W X1 [pi ] W X2 [p 2 ] J [®a]F [a] T£ v [^[pi,p 2 ] + a] (x) . (6) 

The argument = (A,E) denotes collectively the components of the classical fields and their 
canonically conjugate momenta on the initial proper time surface; analytical expressions for these 
are available at T = + [30]. Their temporal evolution is obtained by solving Yang-Mills equations 
with static light front color sources [38]. The W's are the functional density matrices we defined 
previously that obey the JIMWLK equation. To leading log accuracy in x, the Ws for the two 
nuclei can be factorized [39] . The initial spectrum of fluctuations Fo [a] , Gaussian in the quantum 
fluctuations a, has a variance given by the small fluctuation propagator in the Glasma background 
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field as % — > + . In practice, the path integral in a is determined by solving the classical Yang-Mills 
equations repeatedly with the initial conditions at X = + given by 

(7) 



Mnit.+ / d\L K c K a^(x) + c* K a^*(x) 



Here A collectively denotes the quantum fields and their canonical conjugate momenta. The coef- 
ficients c K , with K collectively denoting the quantum numbers labeling the basis of solutions, are 
random Gaussian-distributed complex numbers. Explicit expressions for the small fluctuations and 
their conjugate momenta, denoted here by <%(x) were obtained in ref. [36]. 

The averaging over multiple initializations of the field leads to decoherence of the classical 
fields, coinciding with the development of an equation of state for the (classicaly) conformal field 
theories considered and eventually isotropization and thermalization of the system. In ref. [37], 
based on the formalism outlined here, it was shown explicitly that a longitudinally expanding scalar 
4 theory isotropizes despite a very rapid initial red shift of the longitudinal pressure. This result is 
shown in fig. 2 (left), where the time scale is in lattice units. This gives us hope that the same will 
occur for the gauge theory case as well where numerical work is in progress. 
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Figure 2: Left: Time evolution of the diagonal elements of T^ v and of the trace of the pressure tensor for 
a longitudinally expanding <t> 4 theory. Figure from ref. [37]. Right: Nearside di-hadron yield in high multi- 
plicity collisions from Glasma graphs. The effect of radial flow on this contribution is seen with increasing 
flow velocity /3 (from bottom to top)= 0.0,0.1,0.2,0.25,0.3. Figure from ref. [48]. 

An outstanding problem is how the system evolving with the above initial conditions will 
thermalize. Clearly this has important ramifications for when hydrodynamics is applicable as a 
description of heavy ion collisions. Recently there has been significant progress in understanding 
how thermalization is achieved for a fixed box both from transport computations [40] and classical 
lattice computations [41]. The situation for a longitudinally expanding gauge theory is still fluid 
with a clear resolution likely in the near future [42]. 

An important feature of the Glasma is that it generates long-range correlations that are local- 
ized on scales ~ l/Qs in the transverse plane [43]. Based on the formalism in ref. [39], its feasible 
to study the evolution of these correlations with energy and the the rapidity separation between cor- 
related gluons [44]. Two particle correlations that are azimuthally collimated around Atp ~ and 
are long range in their rapidity separation Atj were observed in high multiplicity p+p collisions by 
the CMS collaboration [45]. Though they appear on the surface similar to phenomena observed in 
A+A collisions at RHIC [46], their origin is qualitatively different. As argued in [47], the effect is 
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due to doubly inclusive 1 /N c 2 suppressed "Glasma" graphs that are enhanced strikingly by a$ (a 
factor of 10 4 - 10 5 for a reasonable range in the coupling) in high multiplicity events for momenta 
< Qs- This enhancement corresponds to effective strong color source densities being probed of 
order ~ 1 / ^fcts in the colliding protons. This explanation was shown to be quantitatively correct in 
ref. [48] from a comparison to the systematics of the per trigger near side yield. It was also shown 
in this work that increasing radial flow changes the behavior of the associated yield as a function 
of trigger particle p T B , in disagreement with data. This result is shown in fig. 2 (right). In contrast, 
the same analysis for the CMS associated yield, in the same kinematic window, for A+A collisions 
shows that the magnitude and shape of the yield with p T s ' agrees with a flow analysis requiring 
a large collective flow velocity. To summarize, the high multiplicity p + p ridge exhibits an az- 
imuthal near side collimation due to a subtle non-trivial manifestation of saturation dynamics in 
multi -particle production. The A+A ridge, while fundamentally also due to the Glasma dynamics 
generating long range correlations, exhibits (a much larger) azimuthal collimation primarily due to 
collective radial flow of hot matter. 
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Figure 3: Left: Initial energy density distribution in three different heavy ion collision events. 
Top:IP-Glasma model. Bottom: MC-KLN model. Figure from ref. [50]. Right: Gluon mul- 
tiplicity (dNg/dy)/ (Npart/2) at T = 0.4fm/c times 2/3 compared to experimental charged particle 
(dN I 'dy) I '(Npart/2) data for y/s = 200GeV Au+Au and y/s = 2.76TeV Pb+Pb collisions as a function of 
A^part- The pale blue and red bands are a collection of the multiplicities for individual events, with the solid 
lines representing the average multiplicity. Figure from ref. [51]. 

As our last topic, we will discuss a phenomenological model of how strongly correlated glu- 
odynamics in nuclear wavefunctions generates, event-by-event, entropy production and flow in 
heavy ion collisions across a range of centralities and energies. This model, called the IP-Glasma 
model, incorporates a number of sources of event-by-event fluctuations in multiparticle produc- 
tion. A significant source of fluctuations are those of the nucleon positions from event-to-event. 
In addition, there are event-by-event fluctuations in the color charge distributions in each of the 
nucleons. Finally, there are fluctuations in the number of gluons produced from event-to-event. 
The color charge fluctuations in the IP-Glasma model are constrained by the IP-Sat model (which 
incorporates energy and impact parameter dependence in the McLerran- Venugopalan model) [49] 
fits to HERA inclusive and diffractive DIS data on e+p collisions, and the limited fixed target e+A 
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DIS data available. Convoluting these nucleon color charge fluctuations with Woods-Saxon dis- 
tributions of nucleons, one constructs Lorentz contracted two dimensional nuclear color charge 
distributions of the incoming nuclei event-by-event. In contrast to prior models, the scale of trans- 
verse event-by-event fluctuations is l/Qs and not the nucleon size. This IP-Glasma model [50, 51] 
employs the fluctuating gluon fields generated by the IP-Sat model to study the event-by-event 
2+1-D Yang-Mills evolution of gluon fields in the Glasma. The corresponding energy density dis- 
tributions in the Glasma vary on the scale l/Qs and are therefore much more highly localized than 
in models that do not include color charge fluctuations. For a comparison of the initial energy den- 
sity distribution in the IP-Glasma model to the CGC motivated MC-KLN model [56] that does not 
include color charge fluctuations and multiplicity fluctuations on size scales 1 / Q$, see fig. 3 (left). 

The IP-Glasma model does an excellent job of describing the single inclusive multiplicity 
as a function of centrality at both RHIC and the LHC. We find that running coupling effects 
are important for a good description of data. (See fig. 3 (right).) We also find that the event- 
by-event IP Glasma does an excellent job of describing the n-particle multiplicity distributions 
measured at RHIC and the LHC. In ref. [52], it was shown in a perturbative computation that 
the n-particle distribution of the "Glasma flux tube" framework of ref. [43] is a negative bino- 
mial distribution (NBD), where the variables are the average multiplicity and the parameter k = 
C(Ac 2 — 1)Q s S±/2ti, with S± the transverse overlap area in the collision and £ a non-perturbative 
parameter to be determined numerically. A Yang-Mills computation of the double inclusive gluon 
distribution [53] shows that this result is non-perturbatively robust with £ ~ 0.3-1.5. The n-particle 
distribution is extracted numerically in ref. [50] from event by event solutions of Yang-Mills equa- 
tions. The distribution obeys the NBD form, with £ ~ 0.2 for smooth distributions at large values of 
Q$S±. These descriptions of the energy and centrality dependence of multiplicity distributions (for 
further examples, see refs. [54]) are strong indications that the CGC provides the right framework 
for entropy production. 

Turning to flow, a sensitive measure is provided by harmonic flow coefficients v„ of the 
anisotropic azimuthal single particle distribution that are related to the initial spatial eccentricities 
£„ by hydrodynamic flow. How efficiently this is done is a measure of the transport properties of 
strongly correlated QCD matter such as the shear and bulk viscosity to entropy density ratios [55] . 
In ref. [5 1], results were presented for £2 through e^. The conclusions that can be drawn are mainly 
that the purely fluctuation driven odd harmonics £3 and £5 from the IP-Glasma model are larger 
than those from the MC-KLN model. In contrast, £2 is smaller than that computed in the MC-KLN 
model, in particular for larger impact parameters. As a consequence, the ratio £2/£3 is smaller than 
in the MC-KLN model, which decreases the ratio of V2/V3 obtained after hydrodynamic evolution, 
thereby making it more compatible with experimental observations. 

Acknowledgments: I thank my collaborators J.-P. Blaizot, A. Dumitru, K. Dusling, T. Epel- 
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hadronic physics. 



8 



Dynamics of strongly correlated gluons at high energies 



Raju Venugopalan 



References 

[1] L.V. Gribov, E.M. Levin, M.G. Ryskin, Phys. Rept. 100, 1 (1983); A.H. Mueller, J-W. Qiu, Nucl. 
Phys. B 268, 427 (1986). 

[2] L.D. McLerran, R. Venugopalan, Phys. Rev. D 49, 2233 (1994); ibid. 49, 3352 (1994); ibid. 50, 2225 
(1994). 

[3] E. Iancu, R. Venugopalan, Quark Gluon Plasma 3, Eds. R.C. Hwa, X.N. Wang, World Scientific, 
hep-ph/0303204; H. Weigert, Prog. Part. Nucl. Phys. 55, 461 (2005); F. Gelis, E. Iancu, J. 
Jalilian-Marian, R. Venugopalan, arXiv: 1002.0333. 

[4] J. Jalilian-Marian, Phys. Rev. D 85, 014037 (2012); G. A. Chirilli, L. Szymanowski and S. Wallon, 
Phys. Rev. D 83, 014020 (2011); T. Altinoluk, A. Kovner, M. Lublinsky and J. Peressutti, JHEP 0903, 
109 (2009); Y. Hatta, E. Iancu, L. McLerran, A. Stasto and D. N. Triantafyllopoulos, Nucl. Phys. A 
764, 423 (2006); F. Gelis and R. Venugopalan, Nucl. Phys. A 779, 177 (2006). 

[5] Jalilian-Marian, A. Kovner, A. Leonidov, H. Weigert, Nucl. Phys. B 504, 415 (1997); ibid., Phys. Rev. 
D 59, 014014 (1999); E. Iancu, A. Leonidov, L.D. McLerran, Nucl. Phys. A 692, 583 (2001); E. 
Ferreiro, E. Iancu, A. Leonidov, L.D. McLerran, Nucl. Phys. A 703, 489 (2002). 

[6] I. Balitsky, Nucl. Phys. B463, 99-160 (1996). 

[7] Y. V. Kovchegov, Phys. Rev. D60, 034008 (1999). 

[8] J. -P. Blaizot, E. Iancu and H. Weigert, Nucl. Phys. A 713, 441 (2003). 

[9] K. Rummukainen and H. Weigert, Nucl. Phys. A 739, 183 (2004). 

[10] A. Dumitru, J. Jalilian-Marian, T. Lappi, B. Schenke and R. Venugopalan, Phys. Lett. B 706, 219 
(2011). 

[11] I. Balitsky and G. A. Chirilli, Phys. Rev. D 77, 014019 (2008); Y. V. Kovchegov and H. Weigert, 
Nucl. Phys. A 784, 188 (2007); G. Beuf, Phys. Rev. D 85, 034039 (2012). 

[12] I. Balitsky and G. A. Chirilli, arXiv: 1207.3844 [hep-ph]. 

[13] E. Iancu and D. N. Triantafyllopoulos, JHEP 1204, 025 (2012). 

[14] J. L. Albacete and Y. V. Kovchegov, Phys. Rev. D 75, 125021 (2007). 

[15] F. D. Aaron et al. [HI and ZEUS Collaboration], JHEP 1001, 109 (2010) [arXiv:091 1.0884 [hep-ex]]. 

[16] J. L. Albacete, J. G. Milhano, P. Quiroga-Arias and J. Rojo, arXiv: 1203. 1043 [hep-ph]. 

[17] J. Jalilian-Marian, Y. V. Kovchegov, Phys. Rev. D70, 114017 (2004). 

[18] F. Dominguez, C. Marquet, B. -W. Xiao, F. Yuan, Phys. Rev. D83, 105005 (2011). 

[19] J. P. Blaizot, F. Gelis, R. Venugopalan, Nucl. Phys. A743, 57-91 (2004). 

[20] A. Dumitru, J. Jalilian-Marian, Phys. Rev. D81, 094015 (2010). 

[21] A. Adare et al. [ PHENIX Collaboration ], Phys. Rev. Lett. 107, 172301 (2011); E. Braidot, 
[arXiv: 1102.0931 [nucl-ex]]. 

[22] C. Marquet, Nucl. Phys. A796, 41-60 (2007). [arXiv:0708.0231 [hep-ph]]. 

[23] J. L. Albacete, C. Marquet, Phys. Rev. Lett. 105, 162301 (2010); K. Tuchin, Phys. Lett. B686, 29-35 
(2010). 



9 



Dynamics of strongly correlated gluons at high energies 



Raju Venugopalan 



[24] F. Dominguez, A. H. Mueller, S. Munier, B. -W. Xiao, [arXiv: 1108. 1752 [hep-ph]]. 
[25] A. Dumitru, J. Jalilian-Marian, Phys. Rev. D82, 074023 (2010). 
[26] A. Stasto, B. -W. Xiao and F. Yuan, arXiv:1109.1817 [hep-ph]. 
[27] K. Kutak and S. Sapeta, arXiv: 1205.5035 [hep-ph]. 

[28] D. Boer, M. Diehl, R. Milner, R. Venugopalan, W. Vogelsang, D. Kaplan, H. Montgomery, S. Vigdor 
et al, [arXiv:l 108.1713 [nucl-th]]; J. L. Abelleira Fernandez [LHeC Study Group Collaboration], 
arXiv: 1206.2913 [physics.acc-ph]. 

[29] M. Froissart, Phys. Rev. 123, 1053 (1961). 

[30] A. Kovner, L. D. McLerran, H. Weigert, Phys. Rev. D52, 3809-3814 (1995); A. Krasnitz, 
R. Venugopalan, Nucl. Phys. B557, 237 (1999). 

[31] T. Lappi, L. McLerran, Nucl. Phys. A772, 200-212 (2006); F. Gelis and R. Venugopalan, Acta Phys. 
Polon. B 37, 3253 (2006). 

[32] S. Mrowczynski, Phys. Lett. B 314, 118 (1993); P. Romatschke and M. Strickland, Phys. Rev. D 68, 
036004 (2003); P. B. Arnold, J. Lenaghan and G. D. Moore, JHEP 0308, 002 (2003). 

[33] P. Romatschke, R. Venugopalan, Phys. Rev. Lett. 96, 062302 (2006); Phys. Rev. D 74, 04501 1 (2006); 
Eur. Phys. J. A29, 71-75 (2006); K. Fukushima, F. Gelis, Nucl. Phys. A 874, 108 (2012). 

[34] S. Y. .Khlebnikov and 1. 1. Tkachev, Phys. Rev. Lett. 77, 219 (1996); L. Kofman, A. D. Linde and 
A. A. Starobinsky, Phys. Rev. D 56, 3258 (1997); R. Micha and 1. 1. Tkachev, Phys. Rev. D 70, 
043538 (2004). 

[35] K. Dusling, T. Epelbaum, F. Gelis and R. Venugopalan, Nucl. Phys. A 850, 69 (2011); T. Epelbaum 
and F. Gelis, Nucl. Phys. A 872, 210 (201 1). 

[36] K. Dusling, F. Gelis and R. Venugopalan, Nucl. Phys. A 872, 161 (201 1). 

[37] K. Dusling, T. Epelbaum, F. Gelis and R. Venugopalan, arXiv: 1206.3336 [hep-ph]. 

[38] A. Krasnitz, R. Venugopalan, Phys. Rev. Lett. 84, 4309 (2000); T. Lappi, Phys. Rev. C67, 054903 
(2003). 

[39] F. Gelis, T. Lappi, R. Venugopalan, Phys. Rev. D78, 054019 (2008); ibid., D78, 054020 (2008); ibid., 
D79, 094017 (2009). 

[40] J. -P. Blaizot, F. Gelis, J. -F. Liao, L. McLerran and R. Venugopalan, Nucl. Phys. A 873, 68 (2012); 
A. Kurkela and G. D. Moore, JHEP 1112, 044 (201 1). 

[41] A. Kurkela and G. D. Moore, arXiv: 1207. 1663 [hep-ph]; S. Schlichting, arXiv: 1207. 1450 [hep-ph]; 
J. Berges, S. Schlichting and D. Sexty, arXiv: 1203 .4646 [hep-ph]. 

[42] J. Berges, J. -P. Blaizot and F. Gelis, arXiv: 1203.2042 [hep-ph]. 

[43] A. Dumitru, F. Gelis, L. McLerran and R. Venugopalan, Nucl. Phys. A 810, 91 (2008). 

[44] K. Dusling, F. Gelis, T. Lappi and R. Venugopalan, Nucl. Phys. A 836, 159 (2010). 

[45] V. Khachatryan et al. [ CMS Collaboration ], JHEP 1009, 091 (2010); W. Li, Mod. Phys. Lett. A 27, 
1230018 (2012). 

[46] J. Adams et al. [STAR Collaboration] Phys. Rev. Lett. 95:152301, (2005); Phys. Rev. C 73, 064907 
(2006); A. Adare et al. [PHENIX Collaboration], arXiv:080 1.4545 [nucl-ex]; B. Alver et al. 
[PHOBOS Collaboration], Phys. Rev. Lett. 104, 062301 (2010). 



10 



Dynamics of strongly correlated gluons at high energies 



Raju Venugopalan 



[47] A. Dumitru, K. Dusling, F. Gelis, J. Jalilian-Marian, T. Lappi, R. Venugopalan, Phys. Lett. B697, 
21-25 (2011). 

[48] K. Dusling and R. Venugopalan, arXiv: 1201.2658 [hep-ph], Phys. Rev. Lett. 108, 262001 (2012). 

[49] J. Battels, K. J. Golec-Biernat and H. Kowalski, Phys. Rev. D 66, 014001 (2002); H. Kowalski and 
D. Teaney, Phys. Rev. D 68, 114005 (2003). 

[50] B. Schenke, P. Tribedy and R. Venugopalan, Phys. Rev. Lett. 108 (2012) 252301, arXiv: 1202.6646 
[nucl-th]. 

[51] B. Schenke, P. Tribedy and R. Venugopalan, arXiv: 1206.6805 [hep-ph]. 
[52] F. Gelis, T. Lappi and L. McLerran, Nucl. Phys. A 828, 149 (2009). 
[53] T. Lappi, S. Srednyak and R. Venugopalan, JHEP 1001, 066 (2010). 

[54] E. Levin and A. H. Rezaeian, Phys. Rev. D 82, 014022 (2010); P. Tribedy and R. Venugopalan, Nucl. 
Phys. A 850, 136 (2011) [Erratum-ibid. A 859, 185 (2011)]; ibid., Phys. Lett. B 710, 125 (2012). 

[55] B. Schenke, J. Phys. G G 38, 124009 (201 1); T. Hirano and Y. Nara, arXiv:1203.4418 [nucl-th]. 

[56] D. Kharzeev, E. Levin and M. Nardi, Phys. Rev. C 71, 054903 (2005); H. -J. Drescher and Y. Nara, 
Phys. Rev. C 76, 041903 (2007). 



11 



